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ABSTRACT 

We present the results of a search for gravitational waves associated with 154 gamma-ray bursts 
(GRBs) that were detected by satellite-based gamma-ray experiments in 2009-2010, during the sixth 
LIGO science run and the second and third Virgo science runs. We perform two distinct searches: a 
modeled search for coalescences of either two neutron stars or a neutron star and black hole; and a 
search for generic, unmodeled gravitational-wave bursts. We find no evidence for gravitational-wave 
counterparts, cither with any individual GRB in this sample or with the population as a whole. For 
all GRBs we place lower bounds on the distance to the progenitor, under the optimistic assumption 
of a gravitational- wave emission energy of lO -2 M c 2 at 150 Hz, with a median limit of 17Mpc. For 
short hard GRBs we place exclusion distances on binary neutron star and neutron star-black hole 
progenitors, using astrophysically motivated priors on the source parameters, with median values of 
16 Mpc and 28 Mpc respectively. These distance limits, while significantly larger than for a search that 
is not aided by GRB satellite observations, are not large enough to expect a coincidence with a GRB. 
However, projecting these exclusions to the sensitivities of Advanced LIGO and Virgo, which should 
begin operation in 2015, we find that the detection of gravitational waves associated with GRBs will 
become quite possible. 

Subject headings: gamma-ray bursts - gravitational waves - compact object mergers 



1. INTRODUCTION 

Gamma-ray bursts (GRBs) are intense flashes of 7-rays 
which are observed approximately once per day and are 
isotropically distributed over the sky (see, e.g. Mcszaros 
2006, and references therein). The variability of the 
bursts on time scales as short as a millisecond indicates 
that the sources are very compact, while the identifica- 
tion of host galaxies and the measurement of redshifts 
for more than 200 bursts have shown that GRBs are of 
extra-galactic origin. 

GRBs are grouped into two broad classes by 
their characteristic duration and spectral hardness 
(Kouvcliotou et al. 1993). Long GRBs (> 2 s, with 
softer spectra), are related to the collapse of massive 
stars with highly rotating cores (see e.g. reviews Modjaz 
2011; Hjorth & Bloom 2011). The extreme core-collapse 
scenarios leading to GRBs result in the formation of a 
stellar-mass black hole with an accretion disk or of a 
highly-magnetized neutron star; for a review see Wooslcy 
(2012) and references therein. In both cases the emission 
of gravitational waves (GWs) is expected, though the 
amount of emission is highly uncertain. 

The progenitors of most short GRBs (< 2 s, with 
harder spectra) are widely thought to be mergers of neu- 
tron star-neutron star or neutron star-black hole binaries 
(see, e.g. Eichlcr et al. 1989; Narayan et al. 1992; Nakar 
2007; Gehrels et al. 2009), though up to a few percent 
may be due to giant flares from a local distribution 
of soft-gamma repeaters (Duncan & Thompson 1992; 
Tanvir et al. 2005; Palmer et al. 2005; Nakar et al. 2006; 
Frcdcriks et al. 2007; Mazets et al. 2008; Chapman et al. 
2009; Hurley et al. 2010). The mergers, referred to here 
as compact binary coalescences, are expected to be 



strong GW radiators (Thorne 1987). The detec- 
tion of gravitational waves associated with a short 
GRB would provide direct evidence that the pro- 
genitor is indeed a compact binary. With such a 
detection it would be possible to measure component 
masses (Finn & Chernoff 1993; Cutler & Flanagan 
1994) and spins (Poisson & Will 1995), constrain 
neutron star equations of state (Vallisneri 2000; 
Flanagan & Hinderer 2008; Hindcrer et al. 2010; 
Read et al. 2009; Lackey et al. 2012; Pannarale et al. 
2011), test general relativity in the strong-field regime 
(Will 2005), and measure calibration- free luminosity dis- 
tances (Schutz 1986; Chernoff & Finn 1993; Dragoljub 
1993; Dalai et al. 2006; Nissanke et al. 2010), which 
allow the measurement of the Hubble expansion and 
dark energy. 

Several searches for gravitational waves associated 
with gamma-ray bursts have been performed using 
data from LIGO and Virgo (Abbott ct al. 2005, 2008b; 
Acernese et al. 2007, 2008). Most recently, data from 
the fifth LIGO science run and the first Virgo science run 
were analyzed to search for coalescence signals or unmod- 
eled gravitational- wave bursts associated with 137 GRBs 
from 2005-2007 (Abbott et al. 2010a,b). No evidence for 
a gravitational-wave signal was found in these searches. 
For GRB 051103 and GRB 070201, short-duration GRBs 
with position error boxes overlapping respectively the 
M81 galaxy at 3.6 Mpc and the Andromeda galaxy (M31) 
at 770 kpc, the non-detection of associated gravitational 
waves ruled out the progenitor object being a compact 
binary coalescence in M81 or M31 with high confidence 
(Abbott et al. 2008a; Abadie et al. 2012b). 

Although it is expected that most GRB progeni- 
tors will be at distances too large for the resulting 
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gravitational- wave signals to be detectable by LIGO and 
Virgo (Bergcr ct al. 2005), it is possible that a few GRBs 
could be located nearby. For example, the smallest ob- 
served redshift to date of an optical GRB afterglow is z = 
0.0085 (~ 36 Mpc) for GRB 980425 (Galama ct al. 1998; 
Kulkarni ct al. 1998; Iwamoto et al. 1998); this would be 
within the LIGO- Virgo detectable range for some pro- 
genitor models. Recent studies (Soderberg et al. 2006; 
Chapman et al. 2007; Le & Dermer 2007; Liang ct al. 
2007; Virgili et al. 2009) indicate the existence of a lo- 
cal population of under-luminous long GRBs with an 
observed rate density approximately 10 3 times that of 
the high-luminosity population. Also, observations sug- 
gest that short-duration GRBs tend to have smaller red- 
shifts than long GRBs (Guetta & Piran 2005; Fox et al. 
2005), and this has led to fairly optimistic estimates 
(Abadie et al. 2010; Leonor et al. 2009) for detecting 
associated gravitational-wave emission. Approximately 
90% of the GRBs in our sample do not have measured 
redshifts, so it is possible that one or more could be much 
closer than the typical ~Gpc distance of GRBs. 

In this paper, we present the results of a search for 
gravitational waves associated with 154 GRBs that were 
detected by satellite-based gamma-ray experiments dur- 
ing the sixth LIGO science run and second and third 
Virgo science runs, which collectively spanned the pe- 
riod from 2009 July 7 to 2010 October 20. We search for 
coalescence signals associated with 26 short GRBs and 
unmodclcd GW bursts associated with 150 GRBs (both 
short and long). The search for unmodeled GW bursts 
targets signals with duration < 1 s and frequencies in the 
most sensitive LIGO/Virgo band, approximately 60 Hz 
— 500 Hz. We find no evidence for a gravitational- wave 
candidate associated with any of the GRBs in this sam- 
ple, and statistical analyses of the GRB sample show no 
sign of a collective signature of weak gravitational waves. 
We place lower bounds on the distance to the progenitor 
for each GRB, and constrain the fraction of the observed 
GRB population at low redshifts. 

The paper is organized as follows. Sec. 2 discusses the 
GW signal models that are used in these searches. Sec. 3 
briefly describes the LIGO and Virgo gravitational-wave 
detectors. Sec. 4 describes the GRB sample during the 
2009-2010 LIGO- Virgo science runs, and Sec. 5 summa- 
rizes the analysis procedure for GW burst signals and for 
coalescence signals. The results are presented in Sec. 6 
and discussed in Sec. 7. We conclude in Sec. 8 with some 
comments on the astrophysical significance of these re- 
sults and the prospects for GRB searches in the era of 
advanced gravitational-wave detectors. 

2. GW SIGNAL MODELS 

As noted above, the progenitors of long GRBs are ex- 
treme cases of stellar collapse, while the most plausible 
progenitors of the majority of short GRBs are mergers 
of a neutron star with either another neutron star or a 
black hole. In this section we review the expected GW 
emission associated with each scenario, and the expected 
delay between the gamma-ray and GW signals. 

2.1. GWs from extreme stellar collapse 

Stellar collapse is notoriously difficult to model. 
It necessitates complex micro-physics and full three- 
dimensional simulations, which take years to complete 



for a single initial state. Many simulations that include 
some, but not all, physical aspects have been performed 
for non-extreme cases of core-collapse supernovae, which 
identified numerous potential GW burst emission chan- 
nels; see Ott (2009) for a review. These models predict 
emission of up to 10~ 8 M c 2 through GWs. Given the 
sensitivity of current GW detectors, such GW emission 
models are not detectable from extra-galactic progeni- 
tors. 

However, in the extreme stellar collapse conditions 
which are necessary to power a GRB, more extreme 
GW emission channels can be considered. Several 
semi-analytical scenarios have been proposed which 
produce up to 10~ 2 M c 2 in GWs, all of which corre- 
spond to some rotational instability developing in the 
GRB central engine (Davies et al. 2002; Fryer et al. 
2002; Kobayashi & Meszaros 2003a; Shibata et al. 
2003; Piro & Pfahl 2007; Corsi & Meszaros 2009; 
Romero et al. 2010). In each model the GWs arc 
emitted by a quadrupolar mass distribution rotating 
around the GRB jet axis. Given the observation of a 
GRB, this axis is roughly pointing at the observer, which 
yields circularly polarized GWs (Kobayashi & Meszaros 
2003b). 

For extreme stellar collapses, the arrival of 7-rays 
can be significantly delayed with respect to the GW 
emission. Delays of up to 100 s can be due to sev- 
eral phenomena: the delayed emission of the relativis- 
tic jet (MacFadyen et al. 2001); sub-luminal propaga- 
tion of the jet to the surface of the star in the collapsar 
model for long GRBs (see for example Aloy ct al. 2000; 
Zhang et al. 2003; Wang & Meszaros 2007; Lazzati ct al. 
2009); and the duration, in the observer's frame, of the 
relativistic propagation of the jet before the onset of the 
prompt 7-ray emission (Vedrenne & Atteia 2009). For 
some GRBs, 7-ray precursors have been observed up 
to several hundred seconds before the main 7-ray emis- 
sion peak (Koshut et al. 1995; Burlon ct al. 2009, 2008; 
Lazzati 2005); the precursor could mark the initial event, 
with the main emission following after a delay. 

2.2. GWs from a compact binary progenitor 

The coalescence of two compact objects is usually 
thought of as a three-step process: an inspiral phase, 
where the orbit of the binary slowly shrinks due to the 
emission of GWs; a merger phase, when the two objects 
plunge together; and a ringdown phase, during which the 
newly created and excited black hole settles into a sta- 
tionary state (Shibata & Taniguchi 2011). As the grav- 
itational waves emitted in the inspiral phase dominate 
the signal-to-noise-ratio in current detectors, we focus 
on that phase only. 1 

We consider compact binaries consisting of two neutron 
stars (NS-NS) or a neutron star with a black hole (NS- 
BH). As the objects spiral together, the neutron star(s) 
are expected to tidally disrupt shortly before they coa- 
lesce, creating a massive torus. The matter in the torus 
can then produce highly relativistic jets, which are sup- 
posedly ejected along the axis of total angular momen- 

1 For high-mass systems the merger and ringdown phases can 
contribute significantly to the signal-to-noisc-ratio with current de- 
tectors. However, given the mass range used in GRB searches, the 
merger and ringdown phases can be ignored. 
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turn. While this picture is supported by recent numerical 
simulations (Foucart ct al. 2011; Rczzolla et al. 2011), it 
has not yet been confirmed by complete simulations, and 
the influence of a tilted BH spin is uncertain. 

Contrary to the long GRB case, the onset of 7-ray 
emission is delayed only up to a few seconds com- 
pared to the GW emission, as there is no dense ma- 
terial retaining the jet and other delay effects are at 
most as long as the GRB duration (Vedrenne & Atteia 
2009). Semi-analytical calculations of the final stages of 
a NS-BH coalescence show that the majority of matter 
plunges onto the BH within ~ls (Davies et al. 2005). 
Numerical simulations of the mass transfer suggest a 
timcscale of milliseconds or a few seconds at maximum 
(Faber et al. 2006; Rosswog 2006; Etienne et al. 2008; 
Shibata & Taniguchi 2008). Therefore, an observer in 
the cone of the collimated outflow is expected to observe 
the gravitational-wave signal up to a few seconds before 
the electromagnetic signal from the prompt emission. 

3. LIGO SCIENCE RUN 6 & VIRGO SCIENCE RUNS 2-3 

The LIGO and Virgo detectors are kilometer-scale, 
power-recycled Michelson interferometers with orthog- 
onal Fabry-Perot arms (Abbott ct al. 2004, 2009a; 
Accadia et al. 2012). They are designed to detect grav- 
itational waves with frequencies ranging from ~ 40 Hz 
to several kHz, with maximum sensitivity near 150 Hz. 
There are two LIGO observatories: one located at Han- 
ford, WA and the other at Livingston, LA. The Hanford 
site houses two interferometers: one with 4 km arms (HI) 
and the other with 2 km arms (H2). The Livingston 
observatory has one 4km interferometer (LI). The two 
observatories are separated by a distance of 3000 km, cor- 
responding to a travel time of 10 ms for light or gravita- 
tional waves. The Virgo detector (VI) is in Cascina near 
Pisa, Italy. The time-of-flight separation between the 
Virgo and Hanford observatories is 27 ms, and between 
Virgo and Livingston is 26 ms. 

A gravitational wave is a spacctime metric perturba- 
tion that is manifested as a time-varying quadrupolar 
strain, with two polarization components. Data from 
each interferometer records the length difference of the 
arms and, when calibrated, measures the strain induced 
by a gravitational wave. These data are in the form of a 
time series, digitized at a sample rate of 16384 Hz (LIGO) 
or 20000Hz (Virgo). 

The sixth LIGO science run was held from 2009 July 
07 to 2010 October 20. During this run, the 4 km HI and 
LI detectors were operated at sensitivities that surpassed 
that of the previous 2005-2007 run, with duty factors of 
52% and 47%. The 2 km H2 detector was not operated in 
2009-2010. The second Virgo science run was held from 
2009 July 07 to 2010 Jan 08 with an improvement in 
sensitivity roughly a factor of 2 over Virgo's first science 
run. The third Virgo science run was held from 2010 
Aug 11 to 2010 Oct 20. The overall Virgo duty cycle 
over these two science runs was 78%. Fig. 1 shows the 
best sensitivities, in terms of noise spectral density, of 
the LIGO and Virgo interferometers during these runs. 
The distance at which the LIGO instruments would ob- 
serve an optimally oriented, optimally located coalescing 
neutron-star binary system with a signal-to-noise-ratio 
of 8 reached about 40 Mpc; for Virgo the same figure of 
merit reached about 20 Mpc. 




frequency [Hz] 



Fig. 1. — Best strain noise spectra from the LIGO and Virgo 
detectors during the 2009-2010 science runs. 

The GEO 600 detector (Grotc ct al. 2008), located 
near Hannover, Germany, was also operational in 2009- 

2010, though with a lower sensitivity than LIGO and 
Virgo. We do not use the GEO data in this search as 
the modest gains in the sensitivity to gravitational-wave 
signals would not have offset the increased complexity 
of the analysis. However, GEO data is used in searches 
for gravitational waves coincident with GRBs occurring 
during periods when only one of the LIGO or Virgo de- 
tectors is operational, such as the period between the 
fifth and sixth LIGO science runs and during summer 

2011. The result of those searches will be reported in a 
future publication. 

4. GRB SAMPLE 

We obtained our sample of GRB triggers from 
the Gamma-ray burst Coordinates Network 2 (GCN) 
(Barthclmy 2008), supplemented by the Swift 3 and 
Fermi 4 trigger pages. This sample of GRB triggers 
came mostly from the Swift satellite (Gehrels et al. 2004) 
and the Fermi satellite (Mcegan et al. 2009), but sev- 
eral triggers also came from other spaceborne experi- 
ments, such as MAXI (Matsuoka et al. 2009), SuperAG- 
ILE (Feroci et al. 2007) and INTEGRAL (Winkler et al. 
2003), as well as from time-of-flight triangulation us- 
ing satellites in the third Interplanetary network (IPN) 
(Hurley ct al. 2009). 

In total there are 404 GRBs in our GRB sample during 
the 2009-2010 LIGO- Virgo science runs. About 10% of 
the GRBs have associated redshift measurements, all of 
them evidently beyond the reach of current GW detec- 
tors. Nevertheless, times around these GRBs have been 
analyzed in case of, for example, a chance association 
with an incorrect host galaxy. 

GRBs that occurred when two or more of the LIGO 
and Virgo detectors were operating in a resonant and 

2 http://gcn.gsfc.nasa.gov/ 

3 http:/ /gen. gsfc.nasa.gov/swift_gnd_ana. html 

4 http:/ /heasarc. gsfc.nasa.gov/W3Browse/fermi/fermigtrig. html 
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stable configuration are analyzed. Data segments which 
arc flagged as being of poor quality are excluded from the 
analysis. In total, 154 GRBs were analyzed, out of which 
150 GRBs were analyzed by the GW burst search, and 
26 short GRBs were analyzed by the coalescence search. 
(As the GW data quality requirements are somewhat dif- 
ferent for the unmodeled burst and coalescence searches, 
4 short GRBs analyzed by the coalescence search could 
not be analyzed by the GW burst search.) 

The classification of GRBs into short and long is some- 
what ambiguous (Gehrels et al. 2006; Bloom et al. 2008; 
Zhang et al. 2009; Horvath et al. 2010). Since binary 
mergers are particularly strong sources of gravitational 
radiation, we make use of a more lenient classification to 
identify GRBs which may originate from a binary merger 
(Zhang et al. 2007, 2009). Our selection is based on the 
Tgo duration (the time interval over which 90% of the to- 
tal background-subtracted photon counts are observed) , 
and on visual inspection of all available lightcurves. 
Specifically, we treat as "short" all GRBs with Tgo < 4 s; 
this choice, rather than the standard 2 s cutoff for short 
GRBs, is to ensure we include those short GRBs in the 
tail of the duration distribution. In addition, some of the 
longer-duration GRBs exhibit a prominent short spike at 
the beginning of the lightcurve and an extended longer 
emission (Norris & Bonncll 2006), suggesting that those 
GRBs might be created by the merger of two compact 
objects. Those GRBs were also treated as short GRBs 
and, where necessary, the trigger time used for the coales- 
cence search was shifted by up to a few seconds to match 
the rising edge of the spike (which should correspond to 
the binary coalescence time). This lenient classification 
ensures a relatively complete sample, at the price of sam- 
ple purity - some of the GRBs we analyze as "short" may 
not have a compact binary progenitor. This impurity is 
acceptable for the purpose of GW detection where we do 
not want to miss a potentially observable GW counter- 
part. The final set of 26 short GRBs is given in Tab. 1. 

A large number of GRBs detected by the IPN are not 
reported by the Gamma-ray burst Coordinates Network; 
the result of a search for GWs associated with those 
GRBs will be reported in a future publication. 

5. SEARCHES FOR GWS ASSOCIATED WITH GRBS 

We perform searches for both unmodeled bursts and 
coalescence signals. We begin this section by describ- 
ing the basic methodology and features common to both 
searches, then briefly present the details of the two anal- 
ysis methods. 

5.1. Search Methodology 

Both search pipelines identify an "on-source" time in 
which to search for an associated GW event. This time 
selection is expected to improve by a factor ~1.5 the 
sensitivity of the search compared to an all-sky / all- 
time search (Kochanek & Piran 1993). For the GW burst 
search, we use the interval from 600 s before each GRB 
trigger to either 60 s or the Tgo time (whichever is larger) 
after the trigger as the window in which to search for a 
GW signal. This conservative window is large enough 
to take into account most plausible time delays between 
a GW signal from a progenitor and the onset of the 
gamma-ray signal, as discussed in Sec. 2.1. This window 
is also safely larger than any uncertainty in the definition 



of the measured GRB trigger time. For cases when less 
early GW data are available, a shorter window starting 
120 s before the GRB trigger time is used. This still cov- 
ers most time-delay scenarios. For the binary coalescence 
search, it is believed that the delay between the merger 
and the emission of 7-rays will be small, as discussed in 
Sec. 2.2. We therefore use an interval of 5 s prior to the 
GRB to 1 s following as the on-source window, which is 
wide enough to allow for uncertainties in the emission 
model and in the arrival time of the electromagnetic sig- 
nal (Abbott et al. 2010b). 

The on-source data are scanned by the search algo- 
rithms to detect possible GW transients (either coa- 
lescence or burst), referred to as "events". For both 
searches the analysis depends on the sky position of the 
GRB. GRBs reported by the Swift satellite have very 
small position uncertainty (^C 1°; see Barthclmy et al. 
2005), and the GW searches need only be performed at 
the reported sky location. For GRBs detected by the 
Gamma-ray Burst Monitor (GBM) on the Fermi satellite 
(Mccgan et al. 2009), however, the sky localization re- 
gion can be large (3> 1°), and detection efficiency would 
be lost if the GW searches only used a single sky loca- 
tion. To resolve this problem, searches for poorly local- 
ized GRBs arc done over a grid of sky positions, covering 
the sky localization region (Was 2011; Was et al. 2012). 
We assume a systematic 68% coverage error circle for the 
Fermi/GBM sky localizations with a radius of 3.2° with 
70% probability and a radius of 9.5° with 30% probabil- 
ity (Connaughton 2011), which is added in quadrature 
to the reported statistical error. 

Each pipeline orders events found in the on-source time 
according to a ranking statistic. To reduce the effect of 
non-stationary background noise, candidate events are 
subjected to checks that "veto" events overlapping in time 
with known instrumental or environmental disturbances 
(Aasi et al. 2012). The surviving event with the highest 
ranking statistic is taken to be the best candidate for 
a gravitational- wave signal for that GRB; it is referred 
to as the loudest event (Brady et al. 2004; Biswas et al. 
2009). To estimate the significance of the loudest event, 
the pipelines also analyze coincident data from a period 
surrounding the on-source data, where we do not expect 
a signal. The proximity of this off-source data to the 
on-source data makes it likely that the estimated back- 
ground will properly reflect the noise properties in the 
on-source segment. The off-source data are processed 
identically to the on-source data; in particular, the same 
data-quality cuts and consistency tests are applied, and 
the same sky positions relative to the GW detector net- 
work are used. If necessary, to increase the background 
distribution statistics, multiple time shifts are applied to 
the data streams from different detector sites, and the 
off-source data re-analyzed for each time shift . 

To determine if a GW is present in the on-source data, 
the loudest on-source event is compared to the distribu- 
tion of loudest off-source events. A p-value is defined 
as the probability of obtaining such an event or louder 
in the onsource, given the background distribution, un- 
der the null hypothesis. The triggers with the smallest 
p-values in the searches are subjected to additional fol- 
lowup studies to determine if the events can be associated 
with some non-GW noise artifact, for example due to an 
environmental disturbance. 
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Regardless of whether a statistically significant signal 
is present, we also set a 90% confidence level lower limit 
on the distance to the GRB progenitor for various signal 
models. This is done by adding simulated GW signals 
to the data and repeating the analyses. These signals, 
which are drawn from astrophysically motivated distri- 
butions described in the following sections, are used to 
calculate the maximum distance for which there is a 90% 
or greater chance that such a signal model, if present in 
the on-source region, would have produced an event with 
larger ranking statistic than the largest value actually 
measured. 

5.2. Search for GW bursts 

The search procedure for GW bursts follows that used 
in the 2005-2007 GRB search (Abbott et al. 2010a). All 
GRBs are treated identically, without regard to redshift 
(if known) , fluence or classification. The on-source data 
are scanned by the X-Pipeline algorithm (Sutton et al. 
2010; Was et al. 2012), which is designed to detect short 
GW bursts, < Is, in the 60 — 500 Hz frequency range. 
X-Pipeline combines data from arbitrary sets of detec- 
tors, taking into account the antenna response and noise 
level of each detector to improve the search sensitivity. 
Time-frequency maps of the combined data streams are 
scanned for clusters of pixels with energy significantly 
higher than that expected from background noise. The 
resulting candidate GW events are characterized by a 
ranking statistic based on energy. We also apply consis- 
tency tests based on the signal correlations measured be- 
tween the detectors, assuming a circularly polarized GW, 
to reduce the number of background events. (The circu- 
lar polarization assumption is motivated by the fact that 
the GRB system rotation axis should be pointing roughly 
at the observer, as discussed in Sec. 2.1.) The stringency 
of these tests is tuned by comparing their effect on back- 
ground events and simulated signal events. The back- 
ground samples are constructed using the ±1.5 hours of 
data around the GRB trigger, excluding the on-source 
time. Approximately 800 time shifts of these off-source 
data are used to obtain a large sample of background 
events. 

To obtain signal samples, simulated signals are added 
to the on-source data. The models of GW emission by 
extreme stellar collapse described in Sec. 2.1 do not pre- 
dict the exact shape of the emitted GW signal. As an 
ad-hoc model, we use the GW emission by a rigidly ro- 
tating quadrupolar mass moment with a Gaussian time 
evolution of its magnitude. For such a source with a 
rotation axis inclined by an angle i with respect to the 
observer the received GW signal is a sine-Gaussian 



h+(t) 
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where the signal frequency fo is equal to twice the rota- 
tion frequency, t is the time relative to the signal peak 
time, Q characterizes the number of cycles for which the 
quadrupolar mass moment is large, -Bgw is the total ra- 
diated energy, and r is the distance to the source. We 



consider two sets of such signals with signal frequencies 
fo of 150 Hz and 300 Hz, which covers the sensitive fre- 
quency band of this GW burst search. The inclination 
angle is distributed uniformly in cos b, with b between 
0° and 5°, which corresponds to the typical jet opening 
angle of - 5° observed for long GRBs (Gal- Yam 2006; 
Racusin et al. 2009). 

Systematic errors are marginalized over in the sensi- 
tivity estimation by "jittering" the simulated signals be- 
fore adding them to the detector noise. This includes 
distributing injections across the sky according to the 
gamma-ray satellites' sky location error box, and jitter- 
ing the signal amplitude, phase, and timing in each de- 
tector according to the given detector calibration errors 
(Accadia et al. 2011; Bartos et al. 2011). This procedure 
also ensures that the consistency tests used in the anal- 
ysis are loose enough to allow for such errors. 

5.3. Search for GWs from a compact binary progenitor 

The core of the coalescence search involves correlating 
the measured data against theoretically predicted wave- 
forms using matched filtering (Hclmstrom 1968). GWs 
from the inspiral phase of a coalescence are modeled by 
post-Newtonian approximants in the band of the de- 
tector's sensitivity for a wide range of binary masses 
(Blanchct 2006). The expected GW signal depends on 
the masses and spins of the NS and its companion (either 
NS or BH), as well as the distance to the source, its sky 
position, its inclination angle, and the polarization angle 
of the orbital axis. Matched filtering is most sensitive 
to the phase evolution of the signal, which depends on 
the binary masses and spins, the time of merger, and a 
fiducial phase. The time and phase can be determined 
analytically. Ignoring spin, we can therefore perform 
matched filtering over a discrete two-dimensional bank 
of templates which span the space of component masses. 
This bank is constructed such that the maximum loss in 
signal-to-noise ratio for a binary with negligible spins is 
3% (Cokelaer 2007; Harry & Fairhurst 2011). For this 
search, as in the 2005-2007 GRB search (Abbott et al. 
2010b), we used "TaylorF2" frequency domain templates, 
generated at 3.5 post-Newtonian order (Blanchet et al. 
1995, 2004). While the spin of the components is ignored 
in the template waveforms, we evaluate the efficiency of 
the search using simulated signals including spin, as de- 
scribed below. 

For each short GRB, the detector data streams are 
combined coherently and searched using the methods de- 
scribed in detail in Harry & Fairhurst (2011). Various 
signal consistency tests are then applied to reject non- 
stationary noise artefacts. These include % 2 tests (Allen 
2005; Hanna 2008), a null stream consistency test, and a 
re-weighting of the signal-to-noise-ratio to take into ac- 
count the values recorded by these tests. This is the first 
coherent search for coalescence signals; it has been found 
to be more sensitive to GW signals than the coincidence 
technique used in previous triggered coalescence searches 
of LIGO and Virgo data (Harry & Fairhurst 2011). Tests 
using the simulations described below have also shown 
that this focused coalescence search is a factor of ^2 

5 X-Pipeline also uses sine-Gaussian signals with fo = 100 Hz 
and non-spinning coalescence signals, as discussed in Sec. 5.3, to 
tune the pipeline. 
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more sensitive to coalescence signals than the unmod- 
eled search described in the previous section, justifying 
the use of a specialized search for this signal type. 

To estimate the efficiency of the search and calculate 
exclusion distances for short GRBs, we draw simulations 
from two sets of astrophysically motivated compact bi- 
nary systems: two neutron stars (NS-NS); and a neu- 
tron star with a black hole (NS-BH). The NS masses are 
chosen from a Gaussian distribution centered at 1.4 M© 
(Kiziltan ct al. 2010; Ozcl et al. 2012) with a width of 
0.2 M Q for the NS-NS case, and a broader spread of 
0.4 M for the NS-BH systems, to account for larger 
uncertainties given the lack of observations for such sys- 
tems. The BH masses are Gaussian distributed with a 
mean of 10 M Q and a width of 6M Q . The BH mass 
is restricted such that the total mass of the system is 
less than 25 Mq. For masses greater than this, the NS 
would be 'swallowed whole' by the BH, no massive torus 
would form, and no GRB would be produced (Duez 2010; 
Ferrari et al. 2010; Shibata & Taniguchi 2011). 

Observed pulsar spin periods and assumptions about 
the spindown rates of neutron stars place the NS spin 
periods at birth in the range of 10 to 140 ms, cor- 
responding to an upper limit on S/m 2 of < 0.04 
(Mandcl & O'Shaughncssy 2010), where S denotes the 
spin of the neutron star and m its mass. However, neu- 
tron stars can be spun up to much higher spins (e.g. 
to 716 Hz (Hessels et al. 2006)), hence we conservatively 
assume a maximum spin of S/m 2 < 0.4 correspond- 
ing to a ~1 ms pulsar. Therefore, the spin magnitudes 
are drawn uniformly from the range [0, 0.4]. For BHs 
the magnitudes are chosen uniformly in the [0, 0.98) 
range (Mandcl & O'Shaughncssy 2010). The spins are 
oriented randomly, with a constraint on the tilt angle 
(the angle between the spin direction of the BH and 
the orbital angular momentum). Since the merger needs 
to power a GRB, a sufficiently massive accretion disk 
around the BH is required. Population synthesis stud- 
ies indicate that the tilt angle is predominantly below 
45° (Bclczynski et al. 2008); numerical simulations show 
that for tilt angles larger than 40° the mass of the disk 
will drop rapidly (Foucart et al. 2011) and BHs with tilt 
angle > 60° will 'swallow' the NS completely, leaving no 
accretion disk to power a GRB (Rantsiou et al. 2008). 
In our simulations we use the weakest of these three con- 
straints and set the tilt angle to be < 60°. 

The outflow from a GRB is most likely to be along the 
direction of the total angular momentum J of the sys- 
tem as discussed in Sec. 2.2. Observations suggest that 
this outflow is confined within a cone, whose half-opening 
angle is estimated to range between several degrees to 
over 60° for short GRBs (see e.g. Burrows ct al. 2006; 
Grupe et al. 2006; Dictz 2011). Under the assumption 
that this cone is centered along the total angular mo- 
mentum J of the system, we chose the inclination angle 
between J and the linc-of-sight to the observer to be dis- 
tributed within cones of half-opening angles 10°, 30°, 45° 
and 90°. The majority of the results quoted in this work 
assume a 30° angle. 

The coalescence time is uniform over the on-source re- 
gion, and the sky position of the GRB is jittered accord- 
ing to the reported uncertainty of the location. 

The quoted exclusion distances are marginalized over 



systematic errors that are inherent in this analysis. First, 
there is some uncertainty in how well our PN templates 
will match real GW signals; we expect a loss in signal- 
to-noise-ratio of up to 10% because of this mis-match 
(Abbott et al. 2009b). Second, there is uncertainty in 
the amplitude calibration of the detectors (Bartos et al. 
2011; Accadia et al. 2011); phase and timing calibration 
uncertainties are also present, but are negligible com- 
pared to other sources of errors. 

An opportunistic search for coalescence signals has also 
been performed on the long GRBs. This search is done to 
conservatively account for uncertainties in the details of 
the short/long GRB classification, and for uncertainties 
in the progenitor model of long GRBs for which an as- 
sociated SN signature was excluded (Gchrels et al. 2006; 
Watson ct al. 2007; Gao ct al. 2010). We use the same 
analysis to check for a coalescence signal associated with 
long GRBs, but do not estimate exclusion distances as 
the compact binary coalescence progenitor model is un- 
likely for long GRBs. 

5.4. Significance of p-value distribution 

In addition to evaluating individual p-valucs, we use 
a weighted binomial test to assess whether the obtained 
set of p-values is compatible with the uniform distribu- 
tion expected from noise only, for both the GW burst 
and coalescence searches. This test looks for deviations 
from the null hypothesis in the 5% tail of lowest p-values 
weighted by the prior probability of detection (estimated 
from the GW search sensitivity). The weighted bino- 
mial test is an extension of the binomial test that has 
been used in previous searches for GW bursts associated 
with GRBs (Abbott ct al. 2008b, 2010a). The combina- 
tion of p-values with prior detection probabilities gives 
more weight to GRBs for which the GW detectors had 
better sensitivity and therefore the detection of a GW 
signal is more likely. The details of this test are given in 
Appendix A. 

The result of the weighted binomial test is a single 
ranking statistic ^weighted- The statistical significance 
of the measured Weighted is assessed by comparing to 
the background distribution of this statistic from Monte 
Carlo simulations with p-values uniformly distributed in 
[0, 1]. This yields the overall background probability of 
the measured set of p-values. 

6. RESULTS 

The coalescence analysis has been applied to search for 
signals in coincidence with 26 short GRBs; the GW burst 
analysis has been applied to 150 GRBs, which include 22 
of the 26 short GRBs analyzed by the coalescence search. 
(As mentioned in Sec. 4, 4 of the short GRBs analyzed 
by the coalescence search could not be analyzed by the 
GW burst search.) The lists of analyzed GRBs classified 
as short and long are given in Tab. 1 and Tab. 2. 

6.1. GW burst search results 

The distribution of p-values for each of the 150 GRBs 
analyzed by the GW burst search is shown in Fig. 2. The 
weighted binomial test yields a background probability 
of 25%. Therefore, the distribution is consistent with no 
GW events being present. 

The smallest p-value, 0.15%, has been obtained for 
GRB 100917A. This GRB was localized on the sky by 
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Fig. 2. — Cumulative p-value distribution from the analysis of 
150 GRBs with the GW burst search. The expected distribution 
under the null hypothesis is indicated by the dashed line. 
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Fig. 3. — Cumulative p-value distribution from the analysis of 26 
short GRBs with the coalescence search. For GRBs where no event 
is observed in the on-source region, we can only place a lower bound 
on the p-value, thus we show two distributions where the upper 
(blue solid line) and lower (green dashed line) bound respectively 
was taken for every GRB. The expected distribution under the null 
hypothesis is indicated by the dashed line. 

Swift, however no redshift measurement is available to 
date. The corresponding GW event was obtained by 
combining data from the HI, LI, and VI detectors. A 
study of the environmental and instrumental channels at 
that time yields potential instrumental causes for this 
event, but is not conclusive. Regardless, the measured 
p-value is not significant as determined by the weighted 
binomial test, so this event is not a candidate for a 
gravitational-wave detection. 



6.2. Coalescence search results 

The distribution of p-values for each of the 26 short 
GRBs analyzed by the coalescence search is shown in 
Fig. 3. The result of the weighted binomial test yields 
a background probability of 8%, corresponding to a 1.8- 
sigma deviation from the null hypothesis. However, as 
we mentioned in section 4, we use a lenient classification 
when deciding if GRBs are treated as short or long for the 
purposes of our analyses. If restrict our short GRB sam- 
ple to the more commonly used criterion, T90 < 2 s, then 
we find a background probability of 3%, corresponding 
to a 2.2-sigma deviation. 

This deviation was due to an event found in coinci- 
dence with GRB 100328A, which produced the smallest 
p-value of 1%, and was the GRB to which the search 
had the second best sensitivity. A followup investiga- 
tion of this candidate determined that it was due to 
a noise artifact in the Hanford instrument, which was 
one of a class of glitches caused by a bad power sup- 
ply which contaminated the length and angular control 
servos. No other noteworthy events were found by this 
search and thus there are no potential gravitational- wave 
candidates. The opportunistic search for coalescence sig- 
nals associated with long GRBs did not yield any candi- 
date that was inconsistent with background noise. 

7. ASTROPHYSICAL INTERPRETATION 

Given that no significant event was found in our anal- 
yses, we place limits on GW emission based on the signal 
models discussed in Sec. 2, and assess the potential of a 
similar search with second-generation gravitational- wave 
detectors. 

7.1. Distance exclusion 

For each GRB we derive a 90% confidence lower limit 
on the GRB progenitor distance for various emission 
models using the methodology described in Sec. 5.1. 

The GW burst search provides lower limits on the 
generic GW burst signal emitted by a rotator described 
in Sec. 5.2 for each GRB. We assume that the source 
emitted -Eqw = 10~ 2 M Q c 2 of energy in gravitational 
waves 6 , that the jet opening angle is 5°, and consider 
emission frequencies of 150 Hz and 300 Hz. The distance 
limits are given in Tab. 1 and Tab. 2, and their histogram 
is shown in Fig. 4. The median exclusion distance is 
D ~ ITMpc^gw/IO^MqC 2 ) 1 / 2 for emission at fre- 
quencies around 150 Hz, where the LIGO- Virgo detector 
network is most sensitive. 

The coalescence search sets lower limits on both the 
NS-NS and NS-BH models described in Sec. 5.3 for each 
short GRB, assuming a jet half-opening angle of 30°. 
The distance limits are given in Tab. 1 and a histogram 
of their values is shown in Fig. 5. The median exclu- 
sion distance for NS-NS (NS-BH) coalescences is 16 Mpc 
(28 Mpc) for the 30° cone. We note that these exclusion 
distances are affected by the choice of signal parameter 
priors in Sec. 5.3; for example, Fig. 6 shows the median 

6 We assume here an astrophysical model of a rotator which 
emits GWs mainly along the rotation axis. In previous searches 
(Abbott et al. 2010a, 2008a) an unphysical isotropic GW emission 
of circularly polarized GWs was used. This change in model in- 
creases the distance exclusions presented here by a factor v/5/2 
relative to previous searches. 
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exclusion distances for half-opening angles of 10°, 30°, 
45°, and 90°. Since the amplitude of a GW signal is 
stronger for a face-on binary, the exclusion distance im- 
proves for smaller half-opening angles. With no restric- 
tion on the opening angle, the 90% exclusion distance 
decreases significantly, as there are orientations which 
will give very little observable GW signal in the detector 
network. 

The GW burst and compact binary coalescence exclu- 
sion distances may be compared to those from all-sky 
searches, which look for GWs without requiring associa- 
tion with a GRB or other external trigger. Figure 7 of 
Abadie et al. (2012a) presents 50%-confidencc exclusion 
energies for the all-sky GW burst search on this same 
data set for an assumed source distance of 10 kpc, with 
a best limit of approximately £gw = 2 x 10 _8 M Q c 2 at 
150 Hz. Rcscaling to our nominal value of 10 -2 M Q c 2 
gives an exclusion distance of ~7 Mpc. Was (2011) per- 
forms a more rigorous comparison that accounts for the 
fraction of events that do not produce GRB triggers due 
to the 7-ray beaming. This indicates that for emission 
opening angles in the 5 — 30° range, the GRB triggered 
search should detect a similar number of GW events com- 
ing from GRB progenitors as that detected by the all-sky 
search - between 0.1 and 6 times the number detected by 
the all-sky search. Furthermore, most of the GW events 
found by the GRB triggered search will be new detections 
not found by the all-sky search, illustrating the value of 
GRB satellites for gravitational-wave detection. 

The NS-NS / NS-BH models used for compact bi- 
nary coalescence exclusions stand on much firmer the- 
oretical ground than the model used for GW burst ex- 
clusions. The amplitude of a coalescence signal is well 
known and depends on the masses and spins of the com- 
pact objects whereas the GW burst energy emitted dur- 
ing a GRB is largely unknown and could be orders of 
magnitude smaller than the chosen nominal value of 
Eqw = 10~ 2 M©c 2 . In the pessimistic scenario that 
GRB progenitors have a comparable GW emission to 
core-collapse supernova, the emitted energy could be as 
low as Eqw ~ 10 _8 M Q c 2 . Such a signal would only 
be observable with current gravitational wave detectors 
from a galactic source. 

7.2. Population exclusion 

As well as a per-GRB distance exclusion, we set an 
exclusion on GRB population parameters by combining 
results from the set of analyzed GRBs. To do this, we use 
a simple population model, where all GRB progenitors 
have the same GW emission (standard sirens), and per- 
form exclusion on cumulative distance distributions. We 
parametrize the distance distribution with two compo- 
nents: a fraction F of GRBs distributed with a constant 
co-moving density rate 7 up to a luminosity distance R, 
and a fraction 1 — F at effectively infinite distance. This 
simple model yields a parameterization of astrophysical 
GRB distance distribution models that predict a uniform 
local rate density and a more complex dependence at red- 
shift > 0.1, as the large redshift part of the distribution 

7 While the distribution of the electromagnctically observed 
GRBs which serve as our triggers needs not be uniform in volume, 
this is a reasonable approximation at the distances to which LIGO- 
Virgo are sensitive. 
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Fig. 4. — Histograms across the sample of GRBs of the distance 
exclusions at the 90% confidence level for circularly polarized sine- 
Gaussian GW burst models at 150 Hz and 300 Hz. We assume an 
optimistic standard siren GW emission of £ GW = 10" 2 M c 2 . Set- 
Tab. 1 and Tab. 2 for the exclusion values for each GRB. 
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Fig. 5. — Histograms across the sample of short GRBs of the 
distance exclusions at the 90% confidence level for NS-NS and 
NS-BH systems. See Tab. 1 for the exclusion values for each short 
GRB. 

is well beyond the sensitivity of current GW detectors. 
The exclusion is then performed in the (F, R) plane. Full 
details of the exclusion method are given in Appendix B. 

The exclusion for GW bursts at 150 Hz with Eqw = 
10~ 2 MqC 2 is shown in Fig. 7, whereas the exclusion for 
the coalescence model for short GRBs is shown in Fig. 8. 
Both exclusions are shown in terms of redshift, where we 
assume a flat ACDM cosmology with Hubble constant 
Hq = 70 kms _1 Mpc _1 , dark matter content FIm = 0.27 
and dark energy content = 0.73 (Komatsu ct al. 
2011). The exclusion at low redshift is dictated by the 
number of analyzed GRBs and at high redshift by the 
typical sensitive range of the search. These exclusions 
assume 100% purity of the GRB sample. For purity p 
the cumulative distribution should be rescaled by 1/p; for 
instance, only one third of our short GRB sample has a 
T90 < 2 s. For comparison, each figure also shows the dis- 
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Fig. 6. — Median exclusion distances of compact binary coales- 
cence sources as a function of half-opening angle, sampled at 10°, 
30°, 45°, and 90°. The medians arc computed over the set of 26 
short GRBs, for both NS-NS and NS-BH, at 90% confidence level. 

tribution of measured GRB rcdshifts, for all Swift GRBs 
(Fig. 7) or for all short GRBs (Fig. 8). While the dis- 
tribution of GRBs with measured rcdshifts includes var- 
ious observational biases compared to the distribution of 
all GRBs detected electromagnetically (and on which we 
perform exclusions), it is clear that the exclusions from 
the current coalescence and GW burst searches are not 
sufficient to put any additional constraint on the nature 
of GRBs. 

While this search for gravitational wave signals in co- 
incidence with observed GRBs was not at the sensitivity 
necessary to detect such coincidences, it is interesting 
to consider the chances of detection with the Advanced 
LIGO/Virgo detectors (Acernese ct al. 2009; Harry ct al. 
2010), which should become operational in 2015. At their 
design sensitivity, these detectors should offer a factor 
of 10 improvement in distance sensitivity to both GW 
burst and coalescence signals, dramatically improving 
the chances to make a gravitational-wave observation of 
an electromagnetically detected GRB. 

In Fig. 7 and Fig. 8 we extrapolate the current exclu- 
sion curves to the advanced detector era, by assuming 
a factor 10 increase in sensitivity of the GW detectors 
and a factor 5 increase in the number of GRBs analyzed 
(equivalent to approximately 2.5 years of live observing 
time at the rate that GRBs are currently being reported) . 
These extrapolations show that detection is quite possi- 
ble in the advanced detector era. Even if a detection is 
not made, targeted gravitational wave searches will allow 
us to place astrophysically relevant constraints on GRB 
population models. 

For long GRBs, the Advanced LIGO/Virgo detectors 
will be able to test optimistic scenarios for GW emission 
- those that produce ~ 10 _2 MqC 2 in the most sensitive 
frequency band of the detectors. The sensitive range for 
these systems will include the local population of sub- 
luminous GRBs that produce the low-redshift excess in 
Fig. 7. We note, however, that GW burst emission with 
significantly lower £qw or a t non-optimal frequencies is 
unlikely to be detectable. 
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Fig. 7. — Cumulative redshift distribution F(R) exclusion from 
the analysis of 150 GRBs with the GW burst search. We exclude 
at 90% confidence level cumulative distance distributions which 
pass through the region above the black solid curve. We assume 
a standard siren sine-Gaussian GW burst at 150 Hz with an en- 
ergy of i?GW = 10 — 2 Mqc 2 . We extrapolate this exclusion to Ad- 
vanced LIGO/Virgo assuming a factor 10 improvement in sensitiv- 
ity and a factor 5 increase in number of GRB triggers analyzed. 
The black dashed curve is the extrapolation assuming the same 
standard siren energy of Eqw = 10~ 2 Mqc 2 and the cyan (gray) 
dashed curve assuming a less optimistic standard siren energy of 
E GW = 10 -4 Mqc 2 (Ott ct al. 2006; Romero et al. 2010). For ref- 
erence, the red staircase curve shows the cumulative distribution of 
measured redshifts for Swift GRBs (.lakobsson ct al. 2006, 2012). 
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Fig. 8. — Cumulative redshift distribution F(R) exclusion from 
the analysis of 26 short GRBs with the coalescence search. Assum- 
ing that all the analyzed short GRBs are NS— BH mergers (NS— NS 
mergers), we exclude at 90% confidence level cumulative distance 
distributions which pass through the region above the black solid 
curve (cyan solid curve). The dashed curves are the extrapolation 
of the solid curves to Advanced LIGO/Virgo, assuming a factor 
10 improvement in sensitivity and a factor 5 increase in number 
of GRB triggers analyzed. For reference, the red staircase curve 
shows the cumulative distribution of measured redshifts for short 
GRBs (Dietz 2011). 
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For short GRBs, a coincident GW detection appears 
quite possible. This conclusion is consistent with simple 
estimates such as that of Metzger & Bcrgcr (2012), who 
estimate a coincident observation rate of 3 yr _1 (0.3 yr~ x ) 
for NS-BH systems (NS-NS systems) with the advanced 
detectors. The precise rate of occurrence will depend on 
the typical masses of the compact objects; we are sen- 
sitive to NS-BH systems at a larger distance than NS- 
NS systems. The distribution of binary component spins 
and the jet opening angle will also affect the received 
GW signal strength. The detection rate will also de- 
pend on the shape of the short GRB cumulative distance 
curve at low redshift. One must also remember that we 
used a very optimistic definition of short GRBs to avoid 
missing a potential signal. It is likely that some of the 
short GRBs that we analyzed for coalescence signals were 
not produced by a compact binary progenitor. Even in 
the case that no coalescence signals are detected in coin- 
cidence with short GRBs in the advanced-detector era, 
it should be possible to place astrophysically interesting 
constraints on the physical characteristics of progenitors 
of short GRBs. 

Finally, we note that these extrapolations carry a num- 
ber of other uncertainties. In particular, the actual per- 
formance of future detectors is unknown. Furthermore, 
the extrapolations depend on how well the sky will be 
covered by gamma-ray satellites in 2015 and later com- 
pared to the present day. 



8. CONCLUSION 

We performed searches for gravitational waves coinci- 
dent with gamma-ray bursts during the 2009-2010 sci- 
ence runs of LIGO and Virgo. In total we analyzed 
154 GRBs using two different analysis methods. A GW 
burst search looked for unmodeled transient signals, as 
expected from massive stellar collapses, and a focused 
search looked for coalescence signals from the merger of 
two compact objects, as expected for short GRBs. We 
did not detect any gravitational wave in coincidence with 
a GRB in cither search. We set lower limits on the dis- 
tance of each GRB with the GW burst search, and of 
the short GRBs with the coalescence search. The median 
exclusion distances are 17M P c(£;gw/10~ 2 M q c 2 ) 1 / 2 at 
150 Hz for the GW burst search and 16 Mpc (28 Mpc) 
for NS-NS (NS-BH) systems for the coalescence search, 
given the priors on the source parameters described in 
Sec. 5. 

These two searches arc more sensitive than the corre- 
sponding all-sky searches of the same data (Abadic ct al. 
2012c,a), due to the more focused analysis possible given 
the trigger time and sky position information provided 
by the GRB satellites. This improvement is as much as 
a factor of ^2 in distance for the GW burst search. Addi- 
tionally, our exclusion distances are greater because each 
source can be presumed to be favorably oriented rela- 



tive to our line of sight, with limits on misalignment set 
by inferences of short and long GRB jet opening angles. 
Further theoretical studies of GRB central engines and 
observational constraints on jet breaks and jet opening 
angles could allow this and future studies to refine their 
constraints a posteriori. Additionally, improved methods 
of classification of GRBs, and in particular of identify- 
ing GRBs with possible binary progenitors with a lower 
false assignment rate, will improve the performance of 
our population estimates. 

The LIGO and Virgo detectors are currently undergo- 
ing a major upgrade, implementing new techniques to 
greatly increase their sensitivity, and are expected to be- 
gin operations by 2015. With these advanced detectors 
our chances to make a coincident GW observation of a 
GRB are good, but depend strongly on the advanced 
detectors running an extended science run at design sen- 
sitivity and the number of GRBs that will be observed 
electromagnetically. Therefore it is of utmost importance 
to have GRB satellites operating during the advanced 
detector era to provide electromagnetic triggers around 
which a more sensitive search for gravitational waves can 
be performed. 
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APPENDIX 
WEIGHTED BINOMIAL TEST 

In a search for GWs associated with GRBs, data corresponding to each GRB are analyzed independently The results 
of these independent analyses need to be combined into a single GW (non-) detection statement, which accounts for 
both the possibility of a single loud GW event or a population of weak GW signals. This weighted binomial test is an 
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extension of the binomial test used to look for an excess of weak gravitational wave signatures in previous searches for 
GW bursts associated with GRBs (Abbott ct al. 2008b, 2010a). 

The binomial test considers the set {pi}i<i<N GRB °f p-values obtained for a population of -/Vgrb analyzed GRBs, 
sorted increasingly. The smallest iVtaii = 0.05./Vgrb of these p-values are used to search for an excess of weak signals. 
The binomial probability, under the null hypothesis, of obtaining at least k events with p-values less than the actual 
fc-th p- value pk is calculated for 1 < k < N ta n and the minimum of these probabilities is used as a detection statistic: 

5 binomial = - log min J2 (TW 1 - P*) N ~ l ■ ( Al ) 

•Sbinomiai looks for a deviation of the p-value distribution when compared to the uniform distribution expected from 
background, in the low p-value region where an excess of weak gravitational wave signals might be observable. However, 
this detection statistic does not take into account the relative a priori GW detection probabilities; that is the sensitive 
volumes of the GW search associated with each GRB trigger, which depends on the GRB sky position and the 
performance of GW detectors at that time. To reduce the contribution of GRBs for which the GW detector sensitivity 
is poor we construct a weighted binomial test (Was 2011) as follows: 

1. Based on the background and sensitivity to simulated signals, compute the distance dk(i) at which the detection 
efficiency is equal to 50% for GRB k and signal emission model i. 

2. Compute the relative volume ratio Rk(i) = dk(i) 3 / max; di(i) 3 for model i compared to the most sensitive GRB. 

3. Average the relative volume ratio over the different models R^ = mean^i?/^). 

4. Sort the penalized p-values Pk/Rk in increasing order, and compute the detection statistic 

^weighted = - log min ( ) TT . (A2) 

For the GW burst search we use the 2 coalescence models and 3 GW burst models given in Sees. 5.2 and 5.3 to 
construct the weighted binomial test, in order to include a range of possible emission models. For the coalescence 
search we use only the 2 coalescence models, which is appropriate for that more focused modeled search. 

POPULATION EXCLUSION METHOD 

A lack of detection can be interpreted individually for each analyzed GRB with an exclusion distance for given 
GW emission models. But the set of analyzed GRBs can also be considered as a whole, to derive constraints on the 
population of GRBs detected by 7-ray satellites. To perform such an exclusion we use a simple population model 
with all GRB progenitors having the same GW emission (standard sirens), and with a distance distribution with two 
components: a fraction F of GRBs distributed with a constant co-moving rate density up to a luminosity distance 
R, and a fraction 1 — F at effectively infinite distance. This simple model yields a parameterization of astrophysical 
GRB distance distribution models that predict a uniform local rate density and a more complex dependence at redshift 
> 0.1, as the large-redshift part of the distribution is well beyond the sensitivity of current GW detectors. 

For this population model we set a frequentist limit on the F and R parameters by excluding all (F, R) which have 
a 90% or greater chance of yielding an event with ranking statistic greater than the largest value actually measured 
for any of the analyzed GRBs. In our computations we assume a flat ACDM cosmology with Hubble constant 
Ho = 70kms _1 Mpc _1 , dark matter content J7m = 0.27 and dark energy content J7a = 0.73 (Komatsu ct al. 2011). 

In practice for each GRB k we measure the efficiency e^(r) as a function of luminosity distance r, for a given 
GW source model of yielding an event with ranking statistic greater than the largest value actually measured. This 
efficiency is integrated over the volume of radius R, where the sources are distributed with constant rate-density. Using 
the volume element for a flat cosmology, 

f- f r . (BD 

dr (i + 2) [(1 + =) 2 + !fVSi M (i + *) s + n A 

we integrate the efficiency as a function of luminosity distance over the considered volume 

E k (R) = J ° T ' d : 1+ * , (B2) 

V ' fR dV_dr_ ' V ; 

JO dr 1 + z 

where the additional l/(l + z) factor accounts for the redshift of the rate. This volume efficiency is the probability for a 
GRB progenitor to yield an event with higher ranking statistic than the value actually measured, under the assumption 
that the GRB has a distance distributed uniformly within the volume of radius R. This can then be extended to a 
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subset of GRBs {ki, . . . , Ivm} all within the local volume of radius R, to construct the probability of at least one of 
them yielding a higher ranking statistic than the measured one: 

E {kl _ kM] (R) = l- ft [l-E k {R]]. (B3) 

ke{k!,...,k M } 

However, our model predicts that a fraction of GRBs 1 — F will originate from distances larger than R, and thus be 
unobservable. For a given fraction F, the distribution of the number J of GRBs in the local volume for a sample of N 
GRBs is binomial, and all subsets {ki, . . . , kj} of [1, NJ have equal probability, given that we assume no knowledge 
of which of the GRBs are in the local volume and which are not. The probability of there being exactly J GRBs in 
the local volume is given by the binomial probability, 

P(J\N)= (^)f j (1-F) n ~ j , (B4) 

and thus the probability of having a given subset of GRBs within R is 

p({k ll ...,kj})=F J (l-F) N - J . (B5) 

We can then obtain the probability that we would have observed a gravitational-wave signal with higher ranking 
statistic than the one actually measured for at least one of the GRBs, as a function of F and i?, by summing over the 
probability of all possible configurations. This is given by 

N 

Ef(R) = Y, E F J (l-F) N - J E {kl _ kj] (R), (B6) 

J=0{fc 1 ,...fe J }c[l,iV] 

and parameters {F 7 R) for which Ef{R) > 0.9 are excluded at 90% confidence. That is, we exclude any cumulative 
distance distribution model that passes through an excluded (F, R) point and which is uniform up to that point. 

This framework can also be expanded to include a mixed sample of GRBs, with a fraction p of GRBs following the 
given standard siren model, and a fraction 1 — p without any significant GW emission. In that case the cumulative 
distance distribution of the GRBs following the standard siren model is excluded whenever E p p(R) > 0.9; that is, the 
exclusion curve is scaled by a 1/p factor compared to the pure sample case. 

RESULT TABLES 



TABLE 1 Short GRB sample and search results 



Exclusion (Mpc) 

UTC network & GW burst at 



GRB name 


time 


RA 


Dec 


time window 


150 Hz 


300 Hz 


NS-NS 


NS-BH 


7-ray detector 


090720B* 


17:02:56 


13 h 31 m 59 s 


-54°48' 


L1V1 


7.1 


3.8 


8.6 


16.0 


GBM 


090802A 


05:39:03 


5 h 37 m 19 s 


34° 05' 


H1L1V1 


7.3 


2.6 


6.5 


11.3 


GBM and IPN 


090815C 


23:21:39 


4 h 17 m 57 s 


-65°57' 


H1L1V1* 


29.8 


12.0 


24.6 


44.3 


BAT 


090820Bt 


12:13:16 


21 h 13 m 02 s 


-18°35' 


H1V1 


12.2 


5.8 


15.1 


26.3 


GBM 


090831 A* 


07:36:36 


9 h 40 m 23 s 


50°58' 


HlVlt 


7.2 


2.1 


4.6 


8.9 


GBM 


090927 


10:07:16 (+1) 


22 h 55 m 42 s 


-70° 58' 


H1L1V1 


16.0 


9.0 


19.8 


35.1 


BAT 


091018* 


20:48:19 


2 h 08 m 46 s 


-57° 33' 


H1V1 






5.2 


10.0 


BAT 


091126A 


07:59:24 


5 h 33 m 00 s 


-19°16' 


H1V1 






13.9 


25.3 


GBM 


091127* 


23:25:45 


2 h 26 m 19 s 


-18°57' 


L1V1* 


5.9 


2.3 


3.1 


4.9 


BAT 


091208B* 


09:49:57 


l h 57 m 39 s 


16°53' 


H1V1 






11.4 


20.6 


BAT 


100111A* 


04:12:49 


16 h 28 m 06 s 


15°32' 


H1L1 


18.8 


8.6 


17.7 


30.4 


BAT 


100206A 


13:30:05 


3 h 08 m 40 s 


13° 10' 


H1L1 


21.0 


8.8 


19.1 


34.1 


BAT 


100213A 


22:27:48 


23 h 17 m 30 s 


43° 22' 


H1L1 


22.4 


10.0 


24.5 


46.3 


BAT 


100216A 


10:07:00 


10 h 17 m 03 s 


35°31' 


H1L1 


29.1 


13.0 


22.7 


40.1 


BAT 


100316B 


08:01:36 


10 h 54 m 00 s 


-45° 28' 


H1L1 






2.1 


3.7 


BAT 


100322B 


07:06:18 


5 h 05 m 57 s 


42°41' 


H1L1 


18.5 


7.1 


14.8 


25.4 


BAT 


100325B* 


05:54:43 


13 h 56 m 33 s 


-79°06' 


H1L1 


21.8 


8.7 


19.0 


34.3 


GBM 


100328A 


03:22:44 


10 h 23 m 45 s 


47° 02' 


H1L1 


28.9 


12.4 


30.1 


51.3 


GBM 


100515A* 


11:13:09 


18 h 21 m 52 s 


27°0T 


H1L1 


38.2 


17.1 


37.1 


64.5 


GBM 


100517D* 


03:42:08 


16 h 14 m 21 s 


-10°22' 


H1L1 


3.4 


2.7 


7.7 


12.1 


GBM 


100628A 


08:16:40 


15 h 03 m 46 s 


-31°39' 


H1L1 


20.6 


8.3 


20.7 


36.7 


BAT 


100717446ft 


10:41:47 


20 h 17 m 14 s 


19°32' 


H1L1 


31.3 


13.2 


26.5 


46.1 


GBM 


100816A 


00:37:51 


23 h 26 m 57 s 


26° 34' 


L1V1 


9.5 


5.8 


6.6 


11.5 


BAT 


100905A 


15:08:14 (-1) 


2 h 06 m 10 s 


14°55' 


H1L1V1° 


17.3 


6.3 


11.5 


19.6 


BAT 


100924A* 


03:58:08 


h 02 m 41 s 


7° 00' 


H1L1V1*° 


29.2 


12.0 


22.8 


39.4 


BAT 


100928A 


02:19:52 (+1) 


14 h 52 m 08 s 


-28°33' 


H1L1V1° 


26.1 


10.1 


20.1 


35.1 


BAT 
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TABLE 1 continued 

Exclusion (Mpc) 

UTC network & GW burst at 

GRB name time RA Dec time window 150 Hz 300 Hz NS-NS NS-BH 7-ray detector 

Information and limits on associated GW emission for each of the analyzed GRBs that were classified by us as short. 
The first four columns are: the GRB name in YYMMDD format or the Fermi/GBM trigger ID for GBM triggers clas- 
sified as a GRB without an available GRB name (see http://heasarc.gsfc.nasa.gov/W3Browse/fermi/fermigbrst.html and 
Pacicsas et al. (2012)); the trigger time (numbers in parentheses denote the time in seconds by which the trigger was shifted 
for the coalescence search following visual inspection of the lightcurve); and the sky position used for the GW search (right 
ascension and declination). Both a " and a * indicate that, although the formal duration of this GRB is longer than 4 s ('), or 
unavailable ("), the GRB was analyzed as a short GRB because of a prominent short spike at the beginning of the lightcurve 
(see Sec. 4). The fifth column gives the gravitational wave detector network used; a * indicates when the shorter on-source 
window starting 120 s before the trigger is used for the GW burst search, and a t when the on-source window is extended to 
cover the GRB duration (T90 > 60s). A indicates the use of only H1L1 data for the burst search, because of data quality 
requirements. Columns 6-9 display the result of the search: the 90% confidence lower limits on the distance to the GRB for 
different waveform models. A standard siren energy emission of i?GW = 10 -2 Mqc 2 is assumed for the circular sine-Gaussian 
GW burst models; these limits are not available for 4 short GRBs which were not analyzed by GW burst search. The last 
column gives the 7-ray detector that provided the sky location used for the search. For GRB 090802, IPN triangulation 
from Konus-WIND, INTEGRAL and Fermi was used to further constrain the sky position. The intersection of the IPN 
and Fermi error regions was used to place search points using the method described in Predoi & Hurley (2012). For this 
GRB, the quoted right ascension and declination corresponds to the centre of the Fermi error region. For iPN localizations 
a complete list of detectors can be found on the project trigger page, http://www.ssl.berkeley.edu/ipn3/masterli.txt. 



TABLE 2 Long GRB sample and search results 



Exclusion (Mpc) 





UTC 






network & 


GW burst at 




GRB name 


time 


RA 


Dec 


time window 


150 Hz 


300 Hz 


7-ray detector 


090709B 


15:07:42 


6 h 14 m 05 s 


64° 05' 


L1V1 


12.4 


6.1 


BAT 


09071 7A 


00:49:32 


5 h 47 m 19 s 


-64° 11' 


HlVl*t 


19.9 


9.7 


GBM 


090719 


01:31:26 


22 h 45 m 04 a 


-67°52' 


H1V1 


10.6 


6.3 


GBM 


090720A 


06:38:08 


13 h 34 m 46 a 


-10°20' 


L1V1 


12.5 


6. 1 


BAT 


090726B 


05:14:07 


16 h 01 m 48 s 


36° 45' 


H1L1V1 


20.2 


6.4 


GBM 


090726 


22:42:27 


16 h 34 m 43 s 


72° 52' 


HlVlt 


17.1 


9.3 


BAT 


090727 


22:42:18 


21 h 03 m 40 a 


64° 56' 


L1V1 1 " 


10.4 


4.9 


BAT 


090727B 


23:32:29 


22 h 53 m 25 a 


-46°42' 


L1V1 


3.3 


1.8 


IPN 


090802B 


15:58:23 


17 h 48 m 04 s 


-71°46' 


H1L1V1 


20.5 


8.3 


GBM 


090807 


15:00:27 


18 h 14 m 57 s 


10° 17' 


HlVlt 


9.8 


5.3 


BAT 


090809 


17:31:14 


21 h 54 m 39 a 


-0°05' 


H1L1V1 


19.2 


6.1 


BAT 


090809B 


23:28:14 


6 h 20 m 60 a 


0°10' 


L1V1 


9.5 


1.9 


GBM 


090810A 


15:49:07 


ll h 15 m 43 s 


-76° 24' 


H1V1 


14.6 


7.1 


GBM 


090814A 


00:52:19 


15 h 58 m 27 s 


25° 35' 


LlVlt 


9.9 


6.1 


BAT 


090814B 


01:21:01 


4 h 19 m 05 s 


60° 35' 


L1V1 


10.9 


5.6 


IBIS 


090814D 


22:47:28 


20 h 30 m 35 s 


45° 43' 


H1L1V1 


17.2 


6.3 


GBM 


090815A 


07:12:12 


2 h 44 m 07 s 


-2°44' 


HlLlVlt 


5.4 


1.3 


GBM 


090815B 


10:30:41 


l h 25 m 40 s 


53° 26' 


H1V1 


10.6 


5.6 


GBM 


090815D 


22:41:46 


16 h 45 m 02 a 


52°56' 


L1V1 


15.7 


6.9 


GBM 


090823B 


03:10:53 


3 h 18 m 07 s 


-17° 35' 


L1V1 


5.9 


2.7 


GBM 


090824A 


22:02:19 


3 h 06 m 35 a 


59°49' 


H1V1 


9.4 


4.8 


GBM 


090826 


01:37:31 


9 h 22 m 28 s 


-0°07' 


H1V1 


2.7 


0,1 


GBM 


090827 


19:06:26 


l h 13 m 44 s 


-50° 54' 


H1V1 


15.1 


8.6 


BAT 


090829B 


16:50:40 


23 h 39 m 57 s 


-9°22' 


HlVlt 


9.0 


4.8 


GBM 


090926B 


21:55:48 


3 h 05 m 14 s 


-38° 60' 


HlLlVlt 


19.1 


7.5 


BAT 


090929A 


04:33:03 


3 h 26 m 47 s 


-7° 20' 


H1V1 


8.7 


1.1 


GBM 


091003 


04:35:45 


16 h 45 m 33 s 


36° 35' 


L1V1 


8.8 


3.2 


LAT 


091017A 


20:40:24 


14 h 03 m ll s 


25°29' 


H1V1 


7.5 


4.9 


GBM 


091018B 


22:58:20 


21 h 27 m 19 s 


-23°05' 


L1V1 


3.6 


1.6 


GBM 


091019A 


18:00:40 


15 h 04 m 07 s 


80° 20' 


H1L1V1 


20.4 


8.2 


GBM 


091020 


21:36:44 


ll h 42 m 54 s 


50° 59' 


L1V1 


9.5 


5.2 


BAT 


091026B 


11:38:48 


9 h 08 m 19 s 


-23°39' 


H1V1 


5.1 


3.0 


GBM 


091030A 


19:52:26 


2 h 46 m 40 s 


21°32' 


HlVlt 


11.5 


1.9 


GBM 


091031 


12:00:28 


4 h 46 m 47 s 


-57° 30' 


HlVlt 


14.2 


6.2 


LAT 


091103A 


21:53:51 


ll h 22 m 24 s 


11°18' 


L1V1 


5.1 


2.2 


GBM 


091109A 


04:57:43 


20 h 37 m 00 s 


-44° 11' 


H1V1 


13.0 


7.3 


BAT 


091109B 


21:49:03 


7 h 31 m 00 s 


-54° 06' 


H1V1 


14.2 


8.7 


BAT 


091115A 


04:14:50 


20 h 31 m 02 s 


71°28' 


H1L1V1* 


14.1 


7.2 


GBM 


091122A 


03:54:20 


7 h 23 m 26 a 


0°34' 


H1V1 


11.9 


5.5 


GBM 


091123B 


01:55:59 


22 h 31 m 16 s 


13°21' 


L1V1* 


9.8 


5,1 


GBM 


091128 


06:50:34 


8 h 30 m 45 s 


1°44' 


H1V1 


6.7 


3.1 


GBM 


091202B 


01:44:06 


17 h 09 m 59 a 


-1°54' 


H1V1 


10.5 


4.6 


GBM 


091202C 


05:15:42 


h 55 m 26 3 


9°05' 


H1V1 


12.6 


6.2 


GBM 


091202 


23:10:04 


9 h 15 m 18 s 


62° 33' 


H1V1 


14.4 


6.5 


IBIS 
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TABLE 2 continued 



Exclusion (Mpc) 











network &£ 


GW burst at 




GRB name 


time 


RA 


Dec 


time window 


150 Hz 


300 Hz 


"-y-ray detector 


091215A 


05-37-26 

UU i U t ■ ij U 


18 h 52 m 59 s 

-LO -> *— UU 


17° 33' 


H1L1* 


18.5 


8.1 


GBM 


091 21 Q A 


11:04:45 


19 h 37 m 57 a 


71°55' 


H1L1V1 


12.4 


6.1 


GBM 


0Q1 220 A 


10-36-50 

-LU . OU . UU 


ll h 07 m 04 s 


4° 49' 


H1L1V1 


21.0 


8.6 


GBM 


0Q1 223R 

U J 1 _ — ■ > L_> 


1 2-1 5-53 

7. _ . _L U . U U 


1 5 h 25 m 04 s 

iU ^U Utt 


54° 44' 


H1L1 


24.3 


9.9 


GBM 


091224A 


08-57-36 


22 h 04 m 40 s 


18° 16' 


H1L1 


16.4 


6. 1 


GBM 


D91 227A 


D7-n3-1 3 

U 1 . UU . -LU 


19 h 47 m 45 a 


2° 36' 


H1L1V1 


19.7 


8.6 


GBM 


i nni m A 


nn-39-49 


20 h 29 m 1fi s 

_ \ ) _ ,7 _L U 


—27° 00' 


H1L1* 


9.1 


1.5 


GBM 


1 nni 03 A 


17:42:32 


7 h 29 m 28 s 


—34° 29' 


H1L1V1 


26.0 


12.1 


IBIS 


1 nni 1 2 A 

± UU -L ± 


10:01:17 


i6 h nn m 33 s 

1U UU UU 


— 75°06' 


H1L1 


15.7 


7.5 


GBM 


100201A 


14:06:17 


8 h 52 m 24 s 


—37° 17' 


H1L1* 


25.9 


10.8 


GBM 


i nn2i 2R 

_L U U £i _L A i_J 


13:11:45 


8 h 57 m 04 s 


32° 13' 


H1L1* 


6.8 


3.7 


GBM 


1 nn2i 

1 U U Jlf lOIJ 


22-58-34 




43° 28' 


H1L1 


15.6 


7.1 


BAT 


100219A 


15:15:46 


10 h 16 m 48 s 


— 12°33' 


H1L1 


17.3 


6.7 


BAT 


1 nn22i A 


08-50-2fi 

UO . 'JU . iiU 


l h 48 m 28 s 


— 17°25' 


H1L1 


29.1 


12.5 


GBM 


i 00225R 


n5-59-n5 

UU . UC7 . UU 


23 h 31 m 24 s 


15° 02' 


H1L1 


20.0 


8.3 


GBM 


1 00225c 


1 3-55-31 


20 h 57 m 04 s 

ZjU U 1 Urt 


0° 13' 


H1L1* 


13.7 


5.8 


GBM 


i 00228R 

_L U U £i £i O i-J 


20:57:47 


7 h 51 m 57 a 


18° 38' 


H1L1 


13.8 


6.8 


GBM 


1 00301 r 


05-21 -4fi 


13 h 27 m 24 s 


19° 50' 


H1L1 


24.0 


9.1 


GBM 


ioosisa 


08-39-1 2 


1 S h 55 m 35 s 

_L U U U UU 


30° 08' 


H1L1 


43.5 


16.9 


GBM 


1 0031 fi A 


02-23-00 


16 h 47 m 48 s 


71°49' 


H1L1 


18.2 


6.9 


BAT 


i 0031 fie 

1UUU 1U V7 


08-57-59 


2 h 09 m 14 s 


—67° 60' 


H1L1 


39.4 


16.5 


BAT 


100324A 


00:21:27 


6 h 34 m 26 s 


—9° 44' 


H1L1 


34.3 


12.6 


BAT 


1 00324R 

1 17 V ' 7 — ^ 1 7 


04-07-3fi 


2 h 38 m 41 s 


—19° 17' 


H1L1 


20.2 


7.3 


IPN 


100325A 


nfi-3fi-ns 


22 h 00 m 57 s 

i_ — UU U 1 


— 26°28' 


H1L1 


36.7 


14.1 


LAT 


i nns2fi A 


r)7-n3-n5 


8 ll 44 m 57 s 


—28° 11' 


H1L1 


13.1 


5.6 


GBM 


ioo°>3iR 


21 -ns-38 

_■ _L i oo ■ oo 


20 h ll m 56 s 


— 11°04' 


H1L1 


14.1 


6.1 


AGILE 


1 nn4m A 

1 V A 7 i V 1 A 


07:07:31 


inhnomi ts 

_L ■ / ' > _L U 


—8° 15' 


HlLlt 


17.6 


6.2 


BAT 


i nn4i n A 


08-31 -57 

UO . O _L . U 1 


8 h 40 m 04 s 


21°29' 


H1L1* 


3.6 


1.5 


GBM 


i nn4i nR 

1 1 7 \) I V 7 I 7 


17;45;4g 


21 h 16 m 59 s 


37° 26' 


H1L1 


28.3 


12.7 


GBM 


1 nn4i s A 

1 V7 V 7 i 1UA 


21 -10-08 

£t±, « _L U i UO 


1 7 h 05 m 25 s 

-L ( UU i(U 


11°27' 


H1L1 


26.5 


11.8 


BAT 


100420B 


00-1 2 0fi 

UU - -L — - uu 


8 h 02 m ll s 


— 5°49' 


H1L1 


32.6 


12.5 


GBM 


i nn42n A 

_L U U rfc a U A 


n5-22-42 


19 h 44 m 21 s 


55°45' 


H1L1 


19.1 


7.5 


BAT 


100423B 


05-51 -25 

UU i U -L ■ iJ u 


7 h 58 m 40 s 


5°47' 


H1L1 


18.6 


6. 1 


GBM 


i 00425 A 


02-50-45 

. UU . rtu 


1 9 h 5fi m 38 s 

LCf UU UO 


— 26°28' 


H1L1 


41.4 


15.6 


BAT 


i nn427A 


08-31 -55 

UO .Jl.uU 


5 h 56 m 41 s 


— 3°28' 


H1L1 


25.5 


11.2 


BAT 


i 00502 A 


08-33-02 

UO . UU . Uii 




18° 23' 


H1L1 


4.0 


2.6 


GBM 


1 00507A 

-LUUUU 1 A 


13:51:15 


h ll m 36 s 


— 79°01' 


H1L1 


23.9 


7.8 


GBM 


mnsnsA 


09:20:42 


5 h 05 m 03 s 

U UU UU 


— 20°45' 


H1L1* 


47.3 


18.2 


BAT 


i 0051 fi A 

1UI/J -LU A 


08-50-41 

UO . UU . _L 


18 h 17 m 38 s 


—8° 12' 


H1L1 


28.9 


11.1 


GBM 


i nn5i fiR 

1 V7 V ' 7 1 UJ—J 


09-30-38 

U t^/ i UU i UO 


19 h 5n m 43 s 

-LU UU i u 


18° 40' 


H1L1 


35.0 


14.3 


GBM 


1 0051 7R 

-LUUU _L 1 i_J 


01 -43-08 

U J. . rtu . UO 


6 h 43 m 43 s 


— 28°59' 


H1L1 


18.6 


7.0 


GBM 


i 0051 7E 

-LUUU _L 1 i—i 


05-49-52 

UU .rtc? . Uii 


gh^^ m 45 s 


4° 26' 


H1L1 


26.8 


10.5 


GBM 


i oo^i 7P 


±u . ± y . oo 


ohorimcrcs 
O OU OtJ 


— 71°52' 


H1L1 


23.9 


10.2 


GBM 


1 0051 7C 

± u u u i- i v_/ 


OS 00 -50 

UU . Uf . UU 


2 h 42 m 31 s 


—44° 19' 


H1L1 


35.9 


13.6 


GBM 


i nn t s2fiR 

_L U U U £i U L_J 


1 9-00-38 

iC7 . UU . UO 


n h n3 m nfi s 

U UU UU 


— 37°55' 


HlLlt 


12.5 


5.4 


BAT 


1 00fi04 A 


UU . UO . O 1 ! 


1 fi h SS m 1 9 s 

1U OJ -LZi 


—73° 11' 


H1L1 


19.3 


9.0 


GBM 


100608A 


0Q-10-0fi 

\J rj . _L U i UU 


2 h 02 m 09 s 


20° 27' 


H1L1 


24.9 


10.4 


GBM 


1 00701 R 


11:45:23 


ohcomofis 
_ ■ >_ _\ ) 


—2° 13' 


H1L1 


14.1 


6.3 


GBM 


1 0070QA 


14:27:32 


9 h 30 m 07 s 


17° 23' 


H1L1 


16.9 


6.5 


GBM 


1 nn71 7372 

1UU i 1 1 iJI i 


08-55-0fi 

UO . UU . UU 


19 h 08 m 14 s 


— 0°40' 


H1L1 


27.3 


10.7 


GBM 


1 nn7i 99S9 


23:44:04 




5° 24' 


H1L1 


15.3 


5.9 


GBM 


100722291 


Ofi-58-24 


2 h 07 m 14 s 


56° 14' 


H1L1* 


17.7 


6.7 


GBM 


1 00725 A 


07-1 2-52 


1 1 h 05 m 52 s 

± _L UU U 


— 26°40' 


HlLlt 


35.5 


14.0 


BAT 


i 00725R 

1UU i Xj ULJ 


11:24:34 


1 9 h 20 m 0fi s 

Li7 ZjU UU 


76° 57' 


HlLlt 


25.9 


10.8 


BAT 


1 nn727A 


05:42:17 


10 h 16 m 44 s 


— 21°25' 


HlLlt 


31.3 


12.6 


BAT 


i 00802 A 


05-45-3fi 

UU . rtU . UU 


h 09 m 55 s 

U Uu UU 


47° 45' 


HlLlt 


36.4 


16.3 


BAT 


1 0080A1 04 


09-9Q-9fi 


1 fihocm cos 
lu uu uZ 


77° 77' 


H1L1 


40.4 


18.0 


GBM 


100814A 


03:50:11 


l h 29 m 54 s 


-17° 59' 


HlLlVlt 


17.3 


7.3 


BAT 


100814351 


08:25:25 


8 h ll m 16 s 


18° 29' 


L1V1* 


14.1 


8.0 


GBM 


100816009 


00:12:41 


6 h 48 m 28 s 


-26°40' 


L1V1* 


6.6 


3.5 


GBM 


100819498 


11:56:35 


18 h 38 m 23 s 


-50°02' 


H1L1V1 


30.1 


12.5 


GBM 


100820373 


08:56:58 


17 h 15 m 09 s 


-18°31' 


H1L1V1 


18.2 


7.2 


GBM 


100823A 


17:25:35 


l h 22 m 49 s 


5°51' 


H1L1V1* 


8.7 


1.0 


BAT 


100825287 


06:53:48 


16 h 53 m 45 a 


-56° 34' 


H1L1V1 


18.7 


7.3 


GBM 


100826A 


22:58:22 


19 h 05 m 43 s 


-32°38' 


L1V1 


4.8 


2.3 


GBM 


100829876 


21:02:08 


6 h 06 m 52 s 


29° 43' 


H1L1V1 


12.3 


4.7 


GBM 
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TABLE 2 continued 



Exclusion (Mpc) 











network & 


GW burst at 




GRB name 


time 


RA 


Dec 


time window 


150 Hz 


300 Hz 


'"y-ray detector 


innQ04A 


01 -33-43 


ll h 31 m 37 s 


— 16° 11' 


L1V1 


13.1 


7.3 


BAT 


i 00905907 

-L VJ yj <3 \J <J <J \J 1 


21 -46-22 


1 7 h 30 m 36 s 

J. 1 'J \J OU 


13° 05' 


H1L1V1 


17.8 


7.5 


GBM 


i nnQfifi A 


13:49:27 


l h 54 m 47 s 


55° 38' 


HlLl*t 


30.5 


12.2 


BAT 


i nfiQi fi A 


18:41:12 


10 h 07 m 50 s 


— 59°23' 


H1L1V1 


8.3 


3.6 


GBM 


100917A 


05 03-25 


19 h 16 m 59 s 


—17° 07' 


HlLlVlt 


18.2 


8.0 


BAT 


1 0091 8863 


20-42-1 8 

^ ' . J__ . ± O 


20 h 33 m 38 s 

ijij O'J UO 


— 45°58' 


H1L1V1 


19.4 


7.8 


GBM 


1 0091 9884 


21:12:16 


10 h 52 m 57 s 


6°01' 


H1L1V1 


19.9 


10.3 


GBM 


1 00922625 


14:59:43 




—25° 11' 


H1V1 


9.2 


4.7 


GBM 




14:17:03 


1 4 hr in m ^Q s 

±■4: JU J3 


— 72°21' 


H1L1 


29.5 


13.1 


GBM 


100926694 


16:39:54 


2 h 54 m 19 s 


-11°06' 


H1L1V1 


12.2 


5.1 


GBM 


100929916 


21:59:45 


12 h 12 m 07 s 


-24° 56' 


H1V1 


10.1 


6.0 


GBM 


101002279 


06:41:26 


21 h 33 m 23 s 


-27° 28' 


H1V1 


9.5 


4.6 


GBM 


101003244 


05:51:08 


ll h 43 m 24 s 


2° 29' 


H1L1V1 


34.1 


13.1 


GBM 


101004426 


10:13:49 


15 h 28 m 52 3 


-43° 59' 


H1L1 


46.8 


17.5 


GBM 


101010190 


04:33:46 


3 h 08 m 45 s 


43° 34' 


L1V1 


9.0 


4.7 


GBM 


101013412 


09:52:42 


19 h 28 m 19 s 


-49° 38' 


H1L1V1 


35.8 


15.0 


GBM 


101015558 


13:24:02 


4 h 52 m 38 s 


15°28' 


H1L1 


29.9 


12.4 


GBM 


101016243 


05:50:16 


8 h 52 m 09 s 


-4° 37' 


L1V1 


9.4 


4.7 


GBM 



Information and limits on associated GW emission for each of the analyzed GRBs that were 
classified as long. The first four columns are: the GRB name in YYMMDD format or the 
Fermi/GBM trigger ID for GBM triggers classified as a GRB without an available GRB name 
(see http://heasarc.gsfc.nasa.gov/W3Browse/fermi/fermigbrst.html and Pacicsas et al. (2012)); 
the trigger time; and the sky position used for the GW search (right ascension and declination). 
The fifth column gives the gravitational wave detector network used; a * indicates when the shorter 
on-source window starting 120 s before the trigger is used, and a t when the on-source window is 
extended to cover the GRB duration (T90 > 60s). Columns 6-7 display the result of the search: 
the 90% confidence lower limits on the distance to the GRB for the circular sine-Gaussian GW 
burst models at 150 Hz and 300 Hz. A standard siren energy emission of -Egw = 10 -2 Mqc 2 is 
assumed. The last column gives the 7-ray detector that provided the sky location used for the 
search. For IPN localizations a complete list of detectors can be found on the project trigger page, 
http: / /www. ssl.berkeley.edu/ipn3/mastcrli.txt. 
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